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ABSTRACT 

Like the Earth, a neutron star (NS) can undergo torque-free precession be- 
cause some piece Aid of its inertia tensor remains tied to the crust's principal 
axes, as opposed to following the crust's angular velocity vector. The (body- 
frame) precession frequency u p is v s AI d /I c , where v s is the NS's spin frequency 
and Ic is the moment of inertia associated with the crustal nuclei, plus any com- 
ponent of the star tightly-coupled to the crust over a timescale less than the spin 
period. For a spinning NS with a relaxed crust, AI d = 6AJq, where AIq is the 
rotational oblateness of a fluid star rotating at spin frequency Q, and b is the 
NS's rigidity parameter. A previous estimate of b by Baym & Pines (1971) gives 
b ~ 10~ 5 for typical NS parameters. Here we calculate the rigidity parameter b, 
and show that it is ~ 40 times smaller than the Baym-Pines estimate. We apply 
this result to PSR 1828-11, an isolated pulsar whose correlated timing residuals 
and pulse shape variations provide strong evidence for precession with a 511-day 
period. We show that this precession period is ~ 250 times shorter than one 
would expect, assuming that: 1) the crust is relaxed (except for the stresses in- 
duced by the precession itself), and, 2) the NS possesses no other source of stress 
that would deform its figure (e.g., a strong magnetic field). We conclude that the 
crust must be under significant stress to explain the precession period of PSR 
B1828-11; such stress arises naturally as the star spins down. Assuming that 
crustal shear stresses do set the precession period, the star's reference angular 
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velocity (roughly, the spin at which the crust is most relaxed) is rs 40 Hz (i.e., 
~ V250 times faster than today's spin), and the weighted average of the crust's 
present spindown strain is a*% e ~ 5 x 10~ 5 . 

We briefly describe the implications of our improved b calculation for other 
precession candidates. 

Subject headings: pulsars: individual(PSR1828-ll) — stars: neutron — stars: ro- 
tation 



1. Introduction 

Stairs et al. (2000) recently reported strong evidence for free precession in the radio 
pulsar PSR 1828-11. This pulsar has period P s = 0.405 s, P s = 6.00 x 10~ 14 , a characteristic 
age of r c =0.11 Myr and an inferred B-field of 5.0 x 10 12 G. Timing residuals for PSR 1828-11 
show strong periodic modulation, with periods 511 and 256 days. A 1009-day periodicity is 
also claimed, but with lower confidence. 

The 511- and 256-day timing residuals are matched by periodic modulations of the 
pulse shape, which strongly supports the interpretation that the modulations are due to 
NS precession, with a wobble amplitude of ~ 3°. In the model of Link & Epstein (2001), 
free precession sets the 511-d timescale, while coupling of the precession to the external 
electromagnetic torque produces variations in the NS's spindown rate that give the main 
contribution to the timing residuals (due to the changing angle between the neutron star's 
magnetic dipole axis and spin direction) and produce a harmonic at 256 d. 

While periodic or quasi-periodic timing residuals have been observed in several other 
isolated pulsars (including B1642-03, the Crab, and Vela) and tentatively interpreted as 
evidence for torque-free precession (see Jones & Anderrson 2001 for a review), PSR 1828-11 
certainly represents the most convincing case for free precession in a neutron star (NS). This 
pulsar has prompted us to re-examine theoretical estimates of a NS's (body-frame) free- 
precession frequency, v v . If the crust is essentially relaxed, the body- frame precession period 
is P p = (ug^Ic/ Ala, where u s is the spin frequency, Iq is the moment of inertia of the crust 
(plus any portion of the NS strongly coupled to the crust on a timescale < and AI d 

is the residual oblateness the NS would have if it were spun down to zero frequency without 
the crust breaking or otherwise relaxing (Pines & Shaham 1972, Munk & MacDonald 1960). 
For a crust that is relaxed at its current spin rate, we can write 



AI d = bAI n , 



(1) 
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where AIn is the piece of NS's inertia tensor due to its spin (i.e., to centrifugal force), 
approximately given by (Jones 2000): 



where Iq is the total moment of inertia for the star when non-rotating. The coefficient b in 
Eq. (1) is sometimes called the NS's "rigidity parameter" (Jones & Andersson 2000). 

This paper is concerned with an accurate calculation of b. For the case of a star with 
constant density and constant shear modulus, the result is b = ^/iV/\E g \, where /i is the 
shear modulus, V is the star's volume, and E g = — |GM 2 /i? is star's gravitational binding 
energy (Love 1944). As an estimate for a neutron star with a liquid core, Baym & Pines 



for a 1.4M NS. The Baym-Pines estimate for b, Eq. (3), has been routinely adopted in the 
NS-precession literature. 

In this paper, we mainly do two things. First, we calculate b for realistic NS structure 
- a solid crust afloat on a liquid core. We solve for the strain field that develops as the NS 
spins down, and find that b is smaller than found by Baym and Pines by a factor of ~ 40. A 
partial explanation of why Eq. (3) yields so large an overestimate is given in the Appendices. 

Second, we use our improved b value to show that PSR 1828-11 is precessing ~ 250 times 
faster than one would expect, assuming its crust is nearly relaxed. Under the hypothesis that 
crust rigidity is the dominant source of the deformation bulge Aid (and that the observed 
modulations are indeed due to precession), we conclude that the crust of PSR 1828-11 is 
not currently relaxed, and that its reference spin z/ s re j (roughly, the spin at which the crust 
is most relaxed) is approximately 40 Hz; i.e., ~ v / 250 times higher than the current spin. 
(In this sense, the precession of PSR B1828-11 is different from the Chandler wobble of the 
Earth, since the Earth is almost relaxed). 

As we will explain in §5, the current strain tensor in the crust <j a b can be written as 
the sum of "spindown" (sd) and "reference" (ref) pieces: cr a b = cr^f + o™ b . The o™ b piece 
is basically the strain that would exist in the crust, even if it were spinning at v s ^ef'i o~ r ^ 
depends on the detailed history of local crust quakes/relaxations, and is thus unknown. The 
spindown strain a^f is the strain induced as the NS spins down from u re f to its current value 
v s . For PSR 1828-11, we show that its weighted-average is a^ c ~ 5 x 10~ 5 . 

Although a a e J is unknown, it seems likely to us that the full strain is larger, on average, 
than the spindown strain. (This is because the full strain a = \ / 0.5cr ab a ab contains contribu- 
tions from all harmonics, while the spindown strain is entirely in the Y 2 q harmonic.) If this 



AI n » 0.3J (z/ s /kHz) 2 , 



(2) 



(1971) took: 




(3) 
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is true, then the NS's crustal breaking strain a max must be larger than the average spindown 
value: a max > ~ 5 x 10~ 5 . Although cr max is poorly constrained, both empirically and 
theoretically, a value this large seems quite reasonable: the usual assumption has been that 
Omax for NS crusts is somewhere in the range 10~ 5 — 10 _1 . (Unlike the crust's shear modu- 
lus, which can be estimated from simple energetics, the breaking strain a max depends on the 
type, density, and propagation of dislocations in the lattice, which are highly uncertain.) 

If significant pinning of superfluid vortex lines to the nuclear lattice or to flux tubes in 
the core occurs, the precession frequency is not set by the material properties of the crust 
(Shaham 1977). In this situation, the pinned superfluid behaves effectively as a large Aid, 
equal to the moment of inertia I p of the pinned superfluid. For example, if the entire crust 
superfluid pins to the solid, I p is ~ 0.01/ ~ Ic- The (body-frame) precession frequency 
should then be of order the spin frequency (or at most ~ 100 times slower, if the crust-core 
coupling is much larger than usually presumed, so Ic ~ Iq). Instead, in PSR 1828-11 and 
the other precession candidates reviewed in §6, the precession frequency is ~ 6 — 8 orders 
of magnitude lower than the spin frequency. Thus claims for long-period NS precession 
appeared at first to conflict with one explanation for large pulsar glitches, wherein superfluid 
vortices do pin to crustal nuclei, and large glitches represent large-scale unpinning events (see, 
e.g., Anderson & Itoh 1975). This conflict was recently resolved by Link & Cutler (2001), 
who showed that, in PSR 1828-11 and other precession candidates, the precessional motion 
itself exerts Magnus forces on the vortices that are ^ 100 times larger than the Magnus forces 
responsible for giant glitches. (The fact that sufficiently large, precession-induced Magnus 
forces would unpin superfluid vortices had been noted previously by Shaham 1986). That 
is, in all likelihood the precessional motion itself immediately unpins the superfluid vortices, 
and keeps them unpinned. Accordingly, in this paper, we set I p = 0. (Of course, pinning 
can persist in precessing stars if the wobble angle is sufficiently small.) 

This paper is concerned with whether crustal stresses alone can give a NS precession 
period of order a year (for spin rates of a few Hz), but we mention that magnetic field stresses 
can also give a precession period of this order if the core is a Type II supeconductor, as shown 
recently by Wasserman (2002), in an extension of work by Mestel and collaborators (Mestel 
& Takhar 1972; Mestel et al. 1981; Nittman & Wood 1981; see all Spitzer 1958). 

The plan of this paper is as follows. In §2 we review the dynamics of NS precession, 
partly to establish notation. In §3 we present our formalism for calculating the rigidity 
parameter b for a NS. In §4 we present the results of this calculation: the value of b and the 
corresponding crustal stresses. Knowing b allows us to estimate the precession frequency of a 
NS with relaxed crust. In §5 we extend this result to a NS with substantially strained crust. 
We apply these results to PSR 1828-11 in §6, where we show that the observed precession 
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frequency is ~ 250 times higher than would be predicted for a relaxed crust. We calculate 
the spindown strain levels that must exist today in the crust of PSR 1828-11, assuming that 
crustal shear stresses are indeed responsible for the observed precession frequency. 

In §7 we briefly discuss the implications of our work for some other NS's that show 
evidence for precession (the Crab, Vela, PSR B1642-03, PSR 2217+47, PSR B0959-54, and 
the possible remnant in SN 1987A). Our conclusions are summarized in §8. 

The Appendices represent our attempts to check and understand our results for b - in 
particular to understand why the Baym-Pines estimate is so much larger than the value we 
obtain. In Appendix A we represent Ala as an integral over crustal shear stresses. The 
result for b we obtain in this way is consistent with our result in §4, and gives more insight 
into how how different stress components and different layers in the crust contribute to the 
final answer. In Appendix B we treat a case where b can be obtained analytically- that of 
a uniform-density, incompressible star with a thin crust of constant shear modulus afloat on 
a liquid core. For this case, we show the actual b is ~ 5 times smaller than the Baym-Pines 
estimate, Eq. (3). In the realistic NS case, b is decreased by additional factor ~ 8, because 
the contributions to AI d from the different stress components and from different layers in 
the crust tend to cancel each other much more nearly than in the case of a thin, uniform 
crust. 



2. Model and Results 

Here we derive the basic equations describing the precession of a NS. In this section, 
and in §§3-6, we will assume that the only sources of non-sphericity in the NS are centrifugal 
force and crustal shear stresses. 

We idealize the precessing NS as having two components: C (for crust) and F (for fluid), 
having angular velocities, and £l F , respectively. and Q F have the same magnitude, f2, 
but different directions. The C piece includes the crustal nuclei, but may also include some 
component of the fluid that is coupled to the crust sufficiently strongly that it effectively 
co-rotates with the crust over timescales shorter than the rotation period 27rf2~ 1 . The F 
piece is that portion of the fluid that is essentially decoupled from the crustal nuclei on this 
timescale; conservation of vorticity then implies that the fluid angular velocity Q F remains 
fixed, and so is aligned with the NS's total angular momentum J a : 

Q a F = QJ a , (4) 

where J a is the unit vector along J a . The angular velocity of the crustal nuclei is Vt^ = ilCl^, 
where fl^ is a unit vector. 
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We write the star's inertia tensor as 

lab = lF,ab + Ic,abi (5) 

where lF,ab an d Ic,ab are the contributions from the fluid interior and the solid crust, re- 
spectively. Following (and somewhat extending) Pines and Shaham (1972), we approximate 
each contribution as the sum of a spherical piece, plus a centrifugal bulge that follows that 
component's angular velocity, plus a deformation bulge (sustained by crustal rigidity) that 
follows some principal crustal axis h a : 



I F ,ab = I Ffi Sab + AI F ^J a J b -^5 ab ) (6) 

+ AI F4 (n a d n b d -^S ab ) 
Ic,ab = Ic,oSab + AI c ,n(n a c n b c -^5 ab ) (7) 
+ AI F , d (n a d n d -±6 ab ). 

where 5 ab is the unit tensor [1, 1, 1]. The sum AI d = AI c> d + AI Fd we call the NS's de- 
formation bulge; it is the non-sphericity the body would retain if it were slowed down to 
zero angular velocity, without the crust breaking or otherwise relaxing. Although AI d is 
ultimately due to crustal shear stresses, the term AI Fd is non-zero (even in the absence of 
pinning), because the crust exerts force on the fluid, both gravitationally and via pressure. 

For a fully relaxed, rotating crust, we expect AI d oc AIq, where AIq = AI f ^q + AIc\n is 
the centrifugal piece of the entire star's inertia tensor, and the coefficient of proportionality 
b = AI d / AIq is the rigidity parameter. For a physical interpretation of b, let us take the 
Earth as an example. The Earth's crust is essentially relaxed. If we could slow the Earth 
down to zero angular velocity without cracking its crust, it would not settle directly into 
a spherical shape, but rather would remain somewhat oblate. This is because the Earth's 
relaxed, zero-strain shape is oblate, and after centrifugal forces are removed, the stresses 
that build up in the crust will tend to push it back toward that relaxed shape. In fact, for 
the Earth, b ~ 0.7 (Munk and MacDonald 1960, p. 40), so stopping the Earth from spinning 
would reduce its oblateness by only ~ 30%. For a NS, as we will see, b ~ 2 x 10 -7 , so halting 
a NS's rotation would decrease its oblateness by a factor of five million. The quantity b 
is so small for a NS because the gravitational energy density far exceeds the crust's shear 
modulus. 

The total angular momentum J a of fluid plus crust is 

Ja = ^(lF,abJ b + Ic,ab&c) (8) 
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Note from Eqs. (7)-(8) that Oq, n a , and J a are all coplanar. Let J a lie along the z-axis, 
and let Oq and h a lie, at some instant, in the x — z plane, as illustrated in Figure 1. The 
angle between the crust's symmetry axis and angular velocity is 6. The angle between the 
angular velocity and the angular momentum is 6. These angles are constants of the motion, 
and satisfy: 



n a = cos 6z a - sin 6 x a 
O a c = cos6z a + sm6x a , 

From Eqs. (9)-(10) we have, up to third-order terms in 9 and 9, 



(9) 
(10) 




r 



9 -oh a . 



(ii) 



From Eq. (11), we see that the precession is a superposition of two motions: n a precesses 
around J a with (inertial frame) precession frequency = Q(l + 9/9), up to terms second 
order in 9 and 9. Likewise, the coefficient of n a in Eq. (11) gives the body-frame precession 
frequency: Q p = 0,(6/6). For O > and AI d > 0, an observer situated above the pole h a and 
fixed in the body frame sees the angular velocity O a circle around h a in the counter-clockwise 
direction. The ratio 6/6 is obtained directly by taking the ^-component of Eq. (8); the result 
is 



6/6 



Ic,o + |A/c,n 



(12) 



3 — l C,d 

where Aid = A/c,d + AIp,d- Since AIc\n and AIc,d are tiny compared to Ic,o, we now ignore 
those terms and drop the subscript "0" from Ic- The (body frame) precession period P p is 
thus 



P p = P s I c /AI d . 



(13) 



Since the pulsar's magnetic dipole moment and emission region are presumably fixed with 
respect to the crust, the angle between the dipole direction and J a varies on the timescale 
P p ] i.e., the body- frame precession period P p is the modulation period for the pulsed radio 
signal. Thus for PSR1828-11, the measured value of AI d /I c is 



Ala 

Ic 



= (v s P p )-> = 9.2 x 10- 9 . 



(14) 



(We have included the factor 511d/P p explicitly in case the precession period is actually 
~ 1009 or 256 days.) 
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3. Elastic Deformation of the Slowing Crust: Formalism 

Our goal in §3-4 is to determine the rigidity parameter b = AI d / AIq, for a NS with a 
relaxed crust. We first show that AI d can be re-expressed as the difference between (i) A/^ 
for a purely fluid NS; and (ii) AI n for the same NS, but now having an elastic crust whose 
relaxed state is spherical. We then describe how we solve for AI n for both cases (i) and (ii), 
using perturbation theory. 

Consider a 2-parameter family of NS's, all with the same mass and composition, and all 
spinning about the z-axis. The two parameters are the NS's actual angular velocity, il, and 
the angular velocity il r at which its crust is relaxed. For convenience, we change variables 
to il 2 and Ail 2 = il 2 — il 2 . The star's inertia tensor I ab is a function of both these variables: 
I a b = labial, Ail 2 ). By axial symmetry, we can write 

I ab = 7 R 2 , Ail 2 )5 ab + I 2 (il 2 r , Ail 2 )(z a z b - ±5 ab ) ■ (15) 

We define the function AI d (il) to be the "residual" I 2 of a NS that is non-rotating, but 
whose crust is relaxed at angular velocity il. Then, by definition, 

Ai d (n) = i 2 (n 2 ,-n 2 ), (16) 

since il 2 = il 2 means that the crust is relaxed at angular velocity il, and Ail 2 = —il 2 means 
that the star is non-rotating. 

Note that J 2 (0, 0) = 0, since il 2 = Ail 2 = means both that the star's crust is relaxed 
and that the star is non-rotating, and hence the star is spherical. Therefore we can write, 
to first order in il 2 and Ail 2 , 



AI d (Q) = I 2 (il 2 ,-il 2 ) (17) 

" n wwr°~ d(AW) lo > (18) 

= I 2 (il 2 ,0)-I 2 (0,il 2 ) . (19) 

This suggests our strategy for computing AI d (il): we will compute both I 2 (il 2 ,0) and 
I 2 (0,il 2 ) and take the difference. Now, I 2 (il 2 ,0) is the I 2 of a star spinning at il, whose 
crust is completely relaxed at that spin. But if the crust is relaxed (i.e., if all shear stresses 
vanish), the shape of the star is the same as for a completely fluid star with the same mass 
and angular velocity. On the other hand, J 2 (0, il 2 ) is the I 2 of a star whose angular velocity 
is il, but whose whose crust is relaxed when spherical (i.e., when it is nonspinning) . 
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Our computational strategy is to start with two nonspinning NS's, (i) and (ii), identical 
except that (i) is treated as a fluid throughout, while (ii) has an elastic crust whose relaxed 
shape is spherical. Now spin them both up to angular velocity Q, calculate I 2 for each, and 
take the difference. The result is AI^ for a NS whose crust is relaxed at spin Q. (More 
precisely, we find Aid(fi), up to fractional corrections of order Q 2 R 3 /GM.) 



3.1. The perturbation equations 

The advantage to the strategy outlined above is that both problems (i) and (ii) can be 
treated with a linear analysis about a spherical background for the stellar structure, with 
centrifugal force acting as a source term that drives the star away from sphericity. The effect 
of centrifugal force is described by the following potential in spherical coordinates (with 
6 = corresponding to the rotation axis): 

S(P C = --Q 2 r 2 (1 - cos 2 9) . (20) 

This potential has both an I = and I — 2, m — piece, 5<p c = <50q + <502, where 

Wo = -^V (21) 

and 

6<j> c 2 (r,9,<j>) = \Y 20 {9,<j>)n 2 r 2 , (22) 

with 

i-il/r- < 23 > 

The I = piece causes spherically-symmetric expansion of the star, while the / = 2 piece gen- 
erates the equatorial bulge. Only the / = 2 piece contributes to I 2 - However, displacements 
due to both I = and I = 2 pieces generate strains in the crust and, hence, are important 
in determining the maximum I 2 spin-change sustainable by the crust without cracking or 
otherwise relaxing. Previous work on spindown-induced crustal strain (e.g., Baym & Pines 
1971; Franco et al. 2000) considered incompressible NS models, where this issue does not 
arise, as the / = piece of the displacement vanishes for incompressible models. We now 
derive the equations governing the response of the star to any 5(p c , and then specialize to 
the I = 2 and I = pieces separately. Note that our treatment will be purely Newtonian. 
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3.1.1. NS response to centrifugal potential 5<p c 

Inside the crust, the perturbed stress tensor (the sum of pressure and shear stresses) 
can be written as 

$r ab = -5p g a b + 2/U a ab , (24) 

where 5p is the pressure perturbation, /i(r) is the crust's shear modulus, and a ab is the crust's 
strain tensor, defined by 

(Tab = \ (V a & + Vbta ~ j^aftV^ ■ (25) 

Here £ a is the displacement of the crust away from its relaxed state. The perturbation 5r ab 
is determined by the condition of hydro-elastic balance, 

V a 5r ab = 5pgr b + pVb (5(f) + 5<j> c 2 ) , (26) 

where f b is the radial unit vector. The perturbation 5(f> of the gravitational potential obeys 

V 2 <50 = AnG5p. (27) 
Finally, the density change arising from the displacement £ a obeys the continuity equation, 

Sp = -V a (p^a) • (28) 

The following bundary conditions arise from Eq. (26): 1) Sr r ± (where "_L" refers to 
components orthogonal to the radial vector f a ) must go to zero at both the crust-core 
boundary and the stellar surface (using the fact that shear stresses vanish just outside the 
crust), and 2) Sr rr = for perturbations that break spherical symmetry (since 5p must vanish 
both in the liquid core-else fluid would flow-and above the star). These conditions can be 
succinctly expressed as 

Sr ab f b = 0, (29) 

at both boundaries. (See Ushomirsky, Cutler & Bildsten 2000 for a more detailed discussion 
of these boundary conditions). 

In the fluid core there are only two first-order equations to solve, for 5(f) and d 8(f)/ dr. 
To derive them, first note that equilibrium between pressure, gravitational, and centrifugal 
forces requires 

V a 5p = -bp V a - pV a [5(f) + 5(f) c ] . (30) 
Projecting Eq. (30) along the horizontal and radial directions, we have 

5p(r) = -p[5(f) + 5(f) c } (31) 
** = -6pg- P ±[6<l> + 6<F\. (32) 
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Taking d/dr of Eq. (31) and subtracting Eq. (32), we have 



so Poisson's equation becomes 



V 2 <50 = AixG{g- ld ^) [6<f> + o0 c ] . (34) 
or 



3.1.2. 1 = 2 deformations 

We now seperate the above equations into spherical harmonics, thereby obtaining a 
system of first-order ODEs. Our treatment follows McDermott et al. (1988); see also 
Ushomirsky, Cutler, & Bildsten (2000, hereafter UCB). Throughout this subsection, it is 
implicit that all scalar perturbations have Y 20 angular dependence; we will generally sup- 
press explicit (/,m) subscripts on perturbed quantities. 

We begin by writing the displacement vector £ a as the sum of radial and tangential 
pieces: 

C = C(r)Y lm f a + ^(O/rVV^™, (35) 

where r a is the unit radial vector and (3 = + 1). Then inside the crust, there are 6 
variables: 

C Ar rr 5t„ d\np 

Zi EE — , z 2 = = Zx— — (36) 

r v v dmr 



o0 R 



2 



c r z 



Ar rr 




V 


V 


Ar rl _ 


Sr r± 


f3p 


Pp 


dz 5 




dlnr 





where 5 and A refer to Eulerian and Lagrangian variations, respectively. The definitions of 
Z1-Z4 and Sr a b are the same as in UCB. The above definition of z 5 , with r 2 in the denominator, 
factors out the dominant behavior of 50/c 2 near r = 0. 

The perturbation equations are (see McDermott et al. 1988; UCB): 

dzi ( a 2 \ 1 a 2 ( s 

1 + 2 — z 1 H z 2 + 6 — z 3 , (37a) 



dlnr \ J 0:3 0:3 

dZ2 = (uV-4V+12r^)z 1+ (v-4^)z 2 (371,) 



dlnr \ a 3 J \ a 3 

+ 6 ( V -6T— ) z 3 + 6z A 

"3/ 
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dz 3 
dlnr 

dz<± 
dliar 



dz 5 
dlnr 

dz 6 
dlnr 



+ —{2z 5 + z 6 ) + 2—\^-, 



1 



-Zi H z 4 , 



V - 6T — ) zi z 2 



a 3 



a 3 



+ — {llaia 2 + lOajj z 3 

a 3 

+ (V-3)z, + —z 5 + —\— 



— Zq, 



-5z 6 + 



UZ* 
c 3 



V m 
- - 4 — | zi 

7 "3 



- ~z 2 + 12^3 

a 3 a 3 



(37c) 
(37d) 



(37e) 
(37f) 

(37g) 



where 



d In q n ~ par dlnv 



dlnr 
<91np 



<91np 



« 2 = r - 



p cilnr 
2/i 



4p 



3p' 



«3 = r + 



7 



fillip 



dlnp 



c 3 



r \ 3 / M 



i?7 Vm, 



3p' 
CM 



i?c 2 ' 



(38) 



In the fluid, we rewrite Eq. (34) in terms of our dimensionless variables z^ and zq: 

dz 5 



dlnr 

dz 6 
dlnr 



UV ( * R 2 Q 2 
-5z 5 z 5 + A — — 

7 V c 



(39) 
(40) 



The boundary conditions are as follows. At the center of the star, d5(p/dr must vanish, 
giving 

z 6 = 0, (41) 
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while Z5 takes on an (unknown) finite value. At both the crust-core boundary (r = r c ) and 
at the top of the crust (r = R), we have, from Eq. (29), 



24 = 0. 

In addition, at the top of the crust (r — R), 

5z 5 + z 6 = 0. (43) 

Finally, at the crust-core boundary 5(f) and d5<f)/dr must be continuous, requiring that z§ 
and z 6 are as well. We thus have 8 equations (two in the core and six in the crust) and 
8 boundary conditions (one at the center, three at the surface, and four at the crust-core 
boundary). For simplicity, we set the outer boundary of the NS at the top of the crust. 
For convenience, we approximate T as 7 in Eqs. (37) (i.e., we do not distinguish d\np/d\n p 
of the background model from the derivative at constant composition). Solutions of these 
equations are described in § 4. 

The I 2 of the NS's inertia tensor is related to the NS's quadrupole moment Q 2 o by 

/ 2 = y^Q 2 o, (44) 

where Q20 is defined by 

Q 20 = J 5p{r) r 4 dr . (45) 
Q20 is related to z 5 at the NS surface by 

5c 2 

Q2o = ~j^z 5 \ r=R . (46) 

As a check on the accuracy of our solutions, we also calculate Q20 directly from the 
following integrals. In the crust, 5p = —V ■ (p£), so 



^ocrust 
V 20 



r R 

2 / p( Zl + 3z 3 )r 4 dr - [rVi]^ (47) 

J r c 

f R U /r\ 5 




(48) 
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In the liquid part of the star, we use Eq. (33) to obtain 



V20 




d\n r. 



(49) 



In a star with both crust and core, we take Q20 — Q20 + Q^o^- m practice, we refine 
our numerical solution until the sum of (48) and (49) agrees with (46) to within required 
tolerance-about one part in 10 10 . 

To find Q20 for a purely fluid star, we solve the same fluid equations, but instead of 
matching at the bottom of the crust, we impose the boundary condition Eq. (43) at r = R. 
Neither solution (i) nor (ii) is really accurate to one part in 10 10 , because we approximate the 
star's outer boundary as being at the top of the crust (ignoring the thin ocean on top), and 
because we approximate T by 7 in Eqs. (37). However since we use the same approximations 
in the solutions to both problems (i) and (ii), the difference is due solely to the presence of 
an elastic crust in model (ii). Therefore we claim we do know difference between (i) and (ii) 
to three significant figures despite the fact that the two quantities agree to seven significant 
figures. This claim is checked in Appendix A, where we obtain an approximate analytical 
expression for b. 



The / = component of the centrifugal potential, 5<pQ, generates a spherically symmetric 
radial displacement field £ a = £ r (r)f a , as well as pressure, density, and gravitational potential 
perturbations. In this subsection, all scalar perturbed quantities are purely functions of r. 

In order to compute £ r for the I = case, we neglect the small shear modulus of the crust 
(effectively taking the entire star to be fluid). Of course, it is the small shear modulus that 
determines AI^, and so in the previous subsection our equations included it. However, our 
present purpose is only to determine (roughly) the / = part of the crustal strain resulting 
from NS spindown, and for this we can neglect p, since it has only a small influence on the 
radial displacement. This approximation is justified by the smallness of the crust's shear 
modulus compared to the pressure. With this approximation, the perturbations in pressure, 
density, and gravitational potential obey 



3.1.3. 1 = deformations 



d5p 



-5pg-p— (6<j> + 6<j>Z) 



(50) 



dr 
Sp 




(51) 



P 



AnGSp. 



(52) 
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Eq. (51) comes from setting ^ = T^. 
We define dimensionless variables 



Vi = —, 2/2 = — (53) 
r p 

5<p dy 3 

2/3 = — , 2/4 



and obtain 



dlnr 



3 h/i 2/2 (54a) 



dlnr \ 7 y 7' 

*> - v(l-l)y 2 -^Ln 2 -l*¥ } ,54b) 



0? In r \ 7/ 1 V r / 3c : 



dlnr 



= y 4 (54c) 



_ UZ± fr \2 

J 2/2-2/4- (54d) 



cilnr C37 

Boundary conditions at r = are obtained from the requirement that the solution be regular 
there. This yields 

4 - ^ yi + ^y 2 = and y 4 = 0. (55) 

At r = R, Ap must vanish and the gravitational potential must match onto its solution in 
empty space (50 oc 1/r), so there 

2/2 - ^2/i = and y 3 + y 4 = , (56) 

while y 3 and 1/4 are continuous at the crust-core interface. 



4. Elastic Deformation of the Slowing Crust: Results 

In this section, we describe the solutions to the elastic perturbation equations derived 
in § 3, and present our results for the value of b = Aid/ In and compare it to the Baym-Pines 
estimate, Eq. (3). 

In order to solve the perturbation equations, we need to specify the background NS 
model. We build NS models using several representative equations of state (EOS), summa- 
rized in Table 1. These models vary somewhat in their treatment of the pressure-density 
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relation in the core and the crust. In particular, model FPS uses the equation of state de- 
scribed in Lorenz, Ravenhall & Pethick (1993); model WS uses EOS UV14+TNI of Wiringa 
et al. (1988) (33) in the core, matched to model FPS in the crust. Models AU and UU use 
equations of state AV14+UVII and UV14+UVII, respectively, matched to EOS of Negele & 
Vautherin (1973) in the crust. In practice, we use tabulations of pressure versus density pro- 
vided with the code RNS 1 . Finally, model nlpoly uses pressure-density relation of an n = 1 
polytrope throughout the star. In all cases, we assume zero temperature, take M = 1.4M , 
and construct the NS model by solving the purely Newtonian equations of hydrostatic bal- 
ance, dp/dr = —pg, and continuity, dM r /dr = Anr 2 p. The resulting NS parameters (radius, 
J , Ic, and AIq/£1 2 ) are shown in Table 1. 

In this paper, we generally assume that Ic equals the moment of inertia of the crustal 
nuclei (i.e., excluding the dripped neutrons in the inner crust), which is roughly 1% of I , 
since estimates of frictional coupling give coupling times between the crust and core well 
in excess of the spin period (see, e.g., Alpar & Sauls 1988). If the stellar magnetic field 
penetrates the core, magnetic stresses will also play a role in the coupling (Abney, Epstein 
& Olinto 1996; Mendell 1998), but the magnetic coupling timescale (estimated as the Alfven 
crossing time for waves in the core) is again much longer than the rotation period for most 
neutron stars. However it is possible that some process does strongly couple the core to crust 
(in which case Ic could be a substantial fraction of Jo), s° we shall indicate how our main 
results scale with Ic/ lo- 
in all models, we place the crust-core boundary at the fiducial density of pj, — 1.5 x 
10 14 g cm -3 (Lorenz, Ravenhall & Pethick 1993). The crust-ocean boundary is placed at 
p = 10 9 g cm~ 3 for all models except for the n = 1 polytrope one, where it is placed at 
10 11 g cm~ 3 in order to ensure numerical stability of our integration of the elastic equations 2 . 

We use the shear modulus p as computed by Strohmayer et al. (1991), via Monte-Carlo 
simulations of both bcc crystals and quenched solids. Their results can be conveniently 
rewritten in terms of the pressure p e of degenerate relativistic electrons, 



where Tcoui = Z 2 e 2 /akT is the ratio of the Coulomb energy of the lattice to the thermal 



1 N. Stergioulas, http://www.gravity.phys.uwm.edu/rns/ 

2 In the n = 1 polytrope model, the density falls of very rapidly with radius in the outer layers-much 
faster than in models with a more realistic EOS. This very rapid variation over many orders of magnitude 
causes numerical problems. However, this truncation docs not affect the final results, as the outer crust 
contains a negligible fraction of the mass. 




(57) 
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energy and Z is the ionic charge. Throughout most of the crust (i.e., except near the top) 
Tcoui > 10 3 , so we for simplicity we ignore the slight dependence of the shear modulus on 
temperature, effectively setting l/T Coul equal to 0. In order to compute the composition 
(A and Z of nuclei, as well as the free neutron fraction) we use the fits of Kaminker et al. 
(1999) to the formalism of Oyamatsu (1993). The resulting run of fi with density is shown 
in Figure 2. For the n = 1 polytrope model we use a fiducial value of fi/p = 10~ 2 . 

With the background model specified as above, we solve the perturbation equations of 
§ 3 using relaxation with adaptive mesh allocation (see, e.g., Press et al. 1992 (24) for a 
description). For the case of I = 2 perturbations in a purely fluid star, as well as for I = 
perturbations, the solution is straightforward. However, for / = 2 perturbations in a star with 
a crust, there are internal boundary conditions at the crust-core interface (Eq. 42), which 
makes application of the standard relaxation method more complicated. In this case, to deal 
with this internal boundary, we transform the independent variable from x — In r to t, such 
that for the fluid part of the star x = (x c — x )t + x$, while in the crust x = x\ — {x\ — x c )t, 
where xq, x±, and x c are the values near the center, at the crust-ocean boundary, and at the 
crust-core boundary, respectively. With this change of the independent variable, both the 
fluid and crust equations can be solved simultaneously on the interval t G [0,1]: the internal 
boundary conditions are eliminated. 

In Figure 3 we show the radial and transverse displacements induced when a relaxed, 
nonrotating star is spun up to Q = c/R. (This is the solution of our perturbation equations 
with (RQ/c) 2 set equal to 1 in our source term Eq. 20.) Of course, Q = c/R is an unreason- 
ably large spin value, but our solutions scale linearly with the source term, i.e., quadratically 
in fl, so in all applications we simply scale the results down to low Q. Recall that, when solv- 
ing the elastic perturbation equations, we imagine taking a spherically symmetric star and 
spinning it up. In this case, the / = part of the centrifugal force acts radially outward, and 
hence the / = component of the radial displacement (lower panel of Figure 3) is positive 
throughout the crust. On the other hand, the / = 2 part of the centrigual force squeezes the 
crust at the poles and pushes it out at the equator. Hence, z\ (top panel of Figure 3), which 
can be thought of as the radial displacement at the poles (times R" 1 ^/ An / 5) , is negative. 
The transverse displacement, z 3 (middle panel Figure 3), changes sign in order to satisfy 
mass continuity. 

Adding the / = and I = 2 perturbations, we can rewrite the strain tensor (Eq. 25) as 
o~ab = (cr + a rr Y 20 )(r a r b - ^e ab ) + a r± f ab + a A (A ab + ^e ab ) . (58) 
Here the tensors e ab , f ab , and A ab are defined by 



= 9ab ~ f a f b , 



(59a) 
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f ab = (3- L r(r a V b Y lm +r b V a Y lm )), (59b) 



^2 



Kb = l Q + 1 ^ aVbYl m + fab, (59c) 



the / = 2 strain components a rr , a r ±, and cta are given by 

2 d^i i 1 2 



Mr) = 3^+3^3 (60a) 

a r± (r) = ^ (60b) 
a A (r) = (3 2 z 3 (60c) 



and the / = piece is 



2 dyi / R1 v 

3 amr 

In Figure 4 we show the radial functions a rr (top panel), <r r ± (second panel), a\ (third panel), 
and cr (bottom panel), for equations of state AU, WS, and nlpoly. Again, these results are 
normalized to (RVt/c) 2 = 1. 

The crust breaks (or deforms plastically) when 

where a max is the yzeW strain. (This is the von Mises criterion; see §6 of UCB for a discussion.) 
From Eq. (58), we get 

2a 2 = - (<7 + a rr y 20 ) 2 + <7? ± -J- (sin 2^) 2 + <7 A ^- sin 4 0, (63) 

2 647T 327T 

where F 20 = (5/16tt) 1 / 2 (3 cos 9 — 1). o" is independent of by symmetry. In Figure 5 we 
plot contours of a on a meridional plane (i.e., plane that slices through the crust of the star 
at a constant longitude 0), and in Figure 6 we show f = 2fia, again on a meridional plane. 
Clearly, the stresses are largest at the base of the crust, where \i is largest, while the strains 
are highest at the top of the crust (where they cost the least, energetically). 

In Table 1 we show the value of b computed as described in § 3 using the solutions to 
the perturbation problem and compare them to the estimates of Baym & Pines (1971) given 
by Eq. (3). Eq. (3) overestimates b by a factor of 15 — 50. To emphasize this difference, and 
also to explore its dependence on p b (the assumed density at the crust-core boundary), we 
show in Figure 7 the ratio of our calculated b to the Baym-Pines estimate as a function of 
Pb. In the two Appendices, we explore why the estimate Eq. (3) is inaccurate. Briefly, when 
a uniform elastic sphere spins down, its strain energy is strongly concentrated toward the 
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center of the star; thus this case (on which Eq. 3 is based) is unsuitable for estimating the 
rigidity of a realistic NS, where all the strain energy is near the surface, in the thin crust. 
Moreover, we show that in a realistic NS crust, contributions to Aid from different stress 
components and different depths cancel each other to a much higher degree than one might 
expect from consideration of a uniform, incompressible crust. 

5. Aid for non-relaxed crusts 

So far in this paper we have calculated the residual oblateness Aid that a spinning NS 
with a relaxed crust retains if it is gently torqued down to zero spin frequency. But a real NS 
that is spinning down (or spinning up, due to accretion) is probably not completely relaxed. 
Indeed, we shall see below that the precession frequency of PSR1828-11 is inconsistent with 
the assumption of a relaxed crust (assuming shear stresses are indeed responsible for this 
NS's AI d ). The crust, because it has finite rigidity, becomes strained as its spins changes. 
(As shown in Fig. 4 and §6, if no cracking or relaxation occurs, then the strain near the 
top of the crust grows to rather large values; e.g., for PSR 1828-11, the strain near the top 
would be ~ 10~ 3 in this case.) In some portions of the crust, the strain could be relaxed in 
localized episodes of crust cracking (starquakes) or plastic flow, while in other regions the 
strain may simply build up. Precisely because the NS crust may have a complex history of 
strain build-up and release, its current strain pattern is impossible to predict exactly. 

Fortunately, however, many of our results can immediately be generalized to this case of 
a non- relaxed crust, by scaling them to the crust's "reference spin" z/ s re j, defined as follows. 
We imagine adjusting the NS's spin frequency v s while keepiing the crust's preferred (zero- 
strain) shape fixed (i.e., without letting the crust crack or otherwise relax). While in general 
there will be no spin at which the crust is completely relaxed, if we plot the resulting 
quadrupole moment Q20 as a function of u s , there would be one spin value, z/ Sjre /, such that 
Q2o(v s ,ref) was identical to the Q20 of a perfect-fluid NS with the same mass and spin, v sre f. 
(Here we are making the natural assumption that the crust's preferred shape is roughly 
oblate, as opposed to prolate.) We define the crustal strain field at v s ,ref to be (J a e J. While 
o~ a e J is in general non-zero, these strains make no net contribution to the NS's quadrupole 
moment. In a rough way, one can think of v s ,ref as the the spin at which the crust is most 
relaxed. As an example, consider a NS that is born with a relaxed crust at spin v s ^ nit and 
never relaxes at all as it spins down to v s , final- Then v s>re f = v s ,init an d o r ^j = 0. If, on the 
other hand, the NS relaxes somewhat (but not completely) as it spins down, then v s ^ef will 
lie somewhere between v s ,init and v s , final, and one would generally expect a a e J 7^ 0. 

The point of these definitions is the following. Imagine a NS with strained crust and spin 
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v s . To determine this NS's precession frequency, we want to know what residual oblateness 
Aid h would have if it were spun down to zero frequency. We can imagine accomplishing 
the spindown in two steps. First, adjust the spin to v sre j. At this spin, the crust has 
some strain <r^ , but o T *j- has no net effect on the NS's oblateness. Second, we spin the 
star down from v s , re f to zero angular velocity. We define a a f to be the extra strain that 
this spindown induces; i.e., a a f = cr a 6 — cr^f, where <r a 5 is the strain in the NS when it is 
nonrotating. Conceptually, we can consider both a r ^j and a a f as independent, small (i.e., 
linear) perturbations, and work to first order in both. To this order, a^f is just the strain 
induced when spinning down a relaxed NS from u s , re f to zero spin, which we solved for in §4. 
Also to this order, the oblateness AI d caused by a ab is the sum of (i) the oblateness due to 
a a f and (ii) the oblateness due to (J^f , but (ii) vanishes, by definition. That is, Aid is just 
the oblateness induced by the spindown part of the strain field, Thus, generalizing Eq. (1), 
we can write 

AI d = bAI n (u s>ref ), (64) 

where AIn(v s%re f) ~ 0.3 Jo(z/ Sjre //kHz) 2 is the oblateness of a fluid star with spin v s ,ref-, and 
b is the same coefficient (i.e., same numerical value) we calculated in § 3. 



6. The precession frequency of PSR 1828-11 

In §2 we showed that the measured precession period of PSR 1828-11 implies that for 
this NS, 

A/ ^l no in-9^ 511c ^ fees 

1 1 8 2 8 = 9.2 x 10 — — . (65) 



1 C \ r V 



For a NS with relaxed crust (i.e., relaxed except for the stresses induced by precessional 
motion itself), the predicted value is 



From the results of §4, we estimate (for our Newtonian NS models): 



A/„//„ * 0.3 ( J!!_V7_«_V , (67) 

' ° VkHz/ \IAM Q ) V 12 kmy ' y ! 



M \ " 2 f R 



Iclh « 0.01 [— ) . (68) 



M \ 3 / R ' 



b w 2 x 10- 7 — — — — . (69) 
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The scaling with NS mass and radius are estimated by taking (i) AIq/I oc MR 2 /(M 2 / R) = 
R 3 /M (i.e., to the ratio of the NS's kinetic and potential energies), (ii) Ic/Io Mq/M oc 
p bot R 2 AR/M oc R 2 Pbot /(gM) oc R 4 /M 2 , and (iii) b oc fi ave V c / (M 2 / R) oc R 2 AR/(M 2 /R) oc 
R 2 (Pbot/ (gPbot)) I (M 2 j R) oc R 5 /M 3 . Here is the mass of the crust, AR is its thickness, 
Vc its volume, pbot and p& ot are the density and pressure at the bottom of the crust, n ave is a 
volume-weighted average of the crustal shear modulus, and g = GM/R 2 . Taking u s = 2.5 Hz, 
we estimate for PSR 1828-11: 

fl^-.Sx^^)^)- 2 ^) 4 . ( 7 0, 

Note that the measured value, Eq. (65), is ~ 250 times larger than our relaxed-crust estimate, 
for our fiducial NS parameters and precession period. And note that if some substantial frac- 
tion of the core is dynamically coupled to the crust (so I c >> 0.01J ), then the discrepancy 
only gets worse. Even if we assume that PSR 1828's precession period is really ~ 1000 days 
and take rather extreme values for the NS mass and radius, M = l.OM and R = 20 km, 
we are still left with a factor ~ 8 discrepancy. We conclude that either our basic picture is 
wrong (i.e., either PSR 1828-11 is not actually precessing, or that some other mechanism 
besides crustal shear stress is responsible for its large Ala) , or that the crust is not relaxed. 
Here we adopt the latter explanation, and pursue its implications. 

The 511-d precession period for PSR B1828-11 implies that the star's reference spin is 
40 Hz. Subsequent spin down to its present spin of 2.5 Hz without significant structural 
relaxation (through quakes or plastic flow) would strain the crust and give the inferred 
AI d /I c = 9.2 x 10 -9 . What are the current strain levels in the crust of 1828-11 in this 
scenario? Compared to 40 Hz, its current spin of 2.5 Hz is very slow. If one continued 
slowing PSR 1828-11 down to zero angular velocity, the spindown-induced crustal stresses 
would increase fractionally by only 0.4%. So (neglecting that 0.4%), the current strains are 
the sum of a) the reference strain a r ^j that this NS would still have if it were spun up to 
v s = 40 Hz and b) the spindown strains cx^f induced by spinning a relaxed NS down from 
v s = 40 Hz to v s = 0. While a r ^ h is practically unknowable, the spindown-induced strains 
a^f are those we found in §4. 

We parametrize the strength of the spindown strain a^f by the scalar quantity a sd , 
defined by 

/I N 1/2 

G = 7>°ab a > ( 71 ) 



.2 

and we define the \i- weighted-average of the crustal spindown strain, o" avc , by 

-sd _ . !v>adV 

avc = r .. ■ i 72 ) 
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We have computed a^ c for our five EOS; the results are in Table 1. They all give roughly 
the same result: a^ e 0.5 for Ril/c = 1. Scaling to the 40 Hz reference angular velocity of 
PSR1828-11 yields 

Again, the total strain a" a (, is the sum of and cx^, where <r^ is determined by the 
detailed evolution of the crust, through structural adjustments, and so without knowing that 
history we can say nothing definitive about a r ^J (except that, by definition, these strains have 
no net effect on the NS's quadrupole moment). However, since a r ^j will contain contributions 
from all harmonics (i.e., all Yi m ), while is wholly I — 2,m — 0, it seems quite likely that 
the NS's average total strain is larger than the average spindown strain. If this is true, the 
average value of spindown portion of the strain places a lower limit on the crustal breaking 
strain: a max ^5x 10 _5 (511d/P p ). In terrestrial solids, strains of order 5 x 10~ 5 are rather 
easy to maintain, so it seems likely that the NS crust is strong enough to sustain the AI d 
implied by the precession frequency. If Ic/Iq is of order unity, the implied average strain 
would be Gave ^ 5x 10~ 3 (511d/P p ), which is still plausibly below the crust's breaking strain. 
We conclude that the neutron star crust is likely strong enough to deform PSR B1828-11 to 
the extent that it will precess with a period of ~ 500 d. 

The spindown strain field suggests that the evolutionary picture of crustal strain might 
be as follows. Suppose the star is born relaxed and rapidly-spinning. As the star spins 
down the strain becomes largest first in the upper crust (see Fig. 5). Cracking begins in 
the upper crust when the yield strain is reached there; however, because relatively little 
strain energy can be stored in this region of small shear modulus, these small quakes would 
not be expect to relax the strain that has developed deeper in the crust. When the yield 
strain is reached deeper, larger (i.e., more energetic) quakes occurs there. Eventually an 
equilibrium is reached in which the strain field is just sub-critical throughout the crust. If 
the Earth is any guide, then numerous small quakes should occur often, while large quakes 
occur rarely, but release most of the accumulated strain energy. These large events could 
excite the precession. We will explore this scenario further in a future publication. 

In this picture, one expects that in the upper crust, a r ^f nearly cancels a^f (so that the 
sum is below the breaking strain), while near the bottom of the crust, where spindown strain 
in lower, the actual strain level a is probably greater than a sd . E.g., if the bottom of the 
crust has not significantly relaxed since the NS's spin was 60 Hz, then the total strain field 
near the bottom would presumably resemble a^f for spindown from v s , re f = 60 Hz. 
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7. Other candidate precessing pulsars 

So far we have focused our discussion on PSR 1828-11; we now apply our results to other 
pulsars showing modulations that may be due to precession. Here we discuss the Crab, Vela, 
PSRs B1642-03, B2217+47, B0959-54, and the possible remnant in SN 1987A. The strength 
of evidence for precession of these pulsars varies from more or less convincing to marginal; 
only one of these candidates (PSR B1642-03) shows strongly periodic and correlated changes 
in both pulse shape and phase similar to those observed in PSR 1828-11. In this section 
we merely summarize all published claims of precession, without attempting to assess their 
validity. Our discussion will be brief, since most these cases have recently been reviewed in 
depth by Jones & Anderson (2001). Note however that our conclusions are often different 
from Jones & Andersson (2001), because our more accurate value of b is nearly two orders 
of magnitude smaller. Our results are summarized in Table 2. 

Before proceeding, we note that two of the precession candidates that we discuss, the 
Crab and Vela pulsars, suffer glitches in spin rate. Glitches are thought to represent angu- 
lar momentum transfer to the crust from the interior superfluid as the array of superfluid 
vortices, usually pinned to nuclei of the crust or magnetic flux tubes in the core, undergoes 
a sudden expansion (see, e.g., Anderson & Itoh 1975; Link & Epstein 1996; Ruderman, Zhu 
& Chen 1998). Even a small amount of vortex pinning is inconsistent with long-period pre- 
cession (Shaham 1977). For the sake of this discussion, we ignore this objection. We refer 
the reader to Link & Cutler (2002) for a discussion of how precession could unpin vortices 
that are intially pinned to the crust. 

7.1. PSR B1642-03 

PSR B1642-03 has v s = 2.5 Hz (the same spin rate as PSR 1828-11), and shows a 
10 3 day periodicity in pulse shape, with modulation amplitude 0.05 (Cordes 1993). This 
period is 120 times shorter than the free-precession period, assuming a relaxed crust. If 
shear stresses are responsible for PSR B1642-03's Aid, then its reference angular velocity is 
~ 27 Hz. 

7.2. PSR B0531+21 (the Crab pulsar) 

The Crab pulsar has P s = 0.0331s; Lyne et al. (1988) observed a phase residual 
with period 20 months, which Jones (1988) interpreted as evidence for precession. Thus 
P s /P p = 6.4 x 10~ 10 , while for our fiducial NS, we predict AI d /I c ~ 5.4 x 10~ 9 This factor 
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~ 8 discrepancy is not large considering the uncertainties in the NS EOS, crust thickness, 
and the mass and radius of the Crab; we conclude that a 20-month precession period is not 
in significant disagreement with our theoretical estimate. This conclusion is contrary to that 
of Jones & Anderson (2001), who concluded that the discrepancy was a factor ~ 700 instead 
of ~ 8. 



7.3. PSR 2217+47 

PSR 2217+47 has P s = 0.538 s, and was originally thought to have a single-component 
pulse profile. However, Suleymanova & Shitov (1994) reported the discovery of a weaker 
second component that varies with a period of ~ 6 — 8 years. Additionally, the braking 
index of 2217+47 has changed significantly between the original observation epoch (1974- 
1984) and the subsequent observations (Shabanova 1990). Suleymanova & Shitov (1994) 
interpret the change in the spindown rate as evidence for free precession with period at least 
comparable to the baseline of observations (i.e., P p ^ 20 years), and attribute the more rapid 
pulse shape variation to a patchy structure of the emission beam. Taking this interpretation 
at face value, the reference angular velocity of the crust is u s , re f <^ 12 Hz (compared to the 
current v s = 1.86 Hz) and a^ c ^ 4 x 10~ 6 . 

7.4. PSR B0959-54 

D'Alessandro & McCulloch (1997) show that the odd moments of the pulse profile of 
B0959-54 are negatively correlated with the variations in the timing residuals. The timing 
residuals appear quasi-periodic, and so they interpret this observation as evidence for free 
precession with P p ^ 2500 days (i.e., greater than or equal to the period of the timing 
residuals). This constraint on P p implies the crust's reference angular velocity is v s , re f ^ 
33 Hz (compared to spin frequency u s = 0.696 Hz), and a^ e ^3x 10~ 5 . 

7.5. PSR B0833-45 (the Vela pulsar) 

The Vela pulsar has P s = 0.089 s; Deshpande & McCulloch (1996) detected intensity 
variations with a 165-d period, which they interpreted as possible evidence of precession. 
Thus P s /P p = 6.2 x 10~ 9 , while for our fiducial NS, we predict AI d /I c « 7 x 10~ 10 . The 
factor ~ 9 discrepancy could be accounted for by Vela having a rather small mass and 
correspondingly large radius (e.g., M « 1.1M and R ~ 18 km), or by Vela having a crust 
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whose reference angular velocity is v s ,ref ~ 32 Hz (i.e., 3 times the current spin). The latter 
would entail an average crustal spindown strain of a^ e ~3x 10~ 5 . 

7.6. Remnant in SN1987A 

Middleditch et al. (2001) have suggested that there is a precessing pulsar in the remnant 
of SN 1987A. They present evidence for the existence of a pulsar (now faded from view) with 
spin period P s = 2.14 ms, and evidence for modulation (which they interpret as precession) 
on a timescale P p ~ 10 3 s, requiring AI^/Ic ~2x 10 -6 . Though failure to observe a pulsar 
in SN1987A since 1993 brings the existence of a pulsar into question, we apply our results 
assuming the pulsar exists. For our fiducial NS parameters, we would expect AI^/Ic ~ 
1.3 x 10 -6 . The factor ~ 1.5 difference is small compared to uncertainties in the NS EOS, 
mass, etc. We conclude (in contrast to Jones & Andersson 2001) that the 10 3 s precession 
timescale is quite reasonable, theoretically. 

7.7. Discussion 

Figure 8 summarizes the observational claims of precession presented in this section. 
In this figure, we plot the observed spin periods of neutron stars against their precession 
claimed periods (or lower limits on P p for PSRs B2217+47 and B0959-54). What limits 
can we place on the precession periods using the interpretation adopted in this paper, i.e., 
that the precession is due to a non-spherical shape of the crust? First, P p = P s (I c /AI d ) = 
P s (Ic/bAI n ) 

^ v s ,ref (f° r §i ven v s)- Dashed lines in Figure 8 show the above relation for 
V^re/ — 1 ms an d 20 ms, respectively. A pulsar lying below the 1 ms line would have a 
reference spin frequency faster than 1 kHz. Note that none of the pulsars with claims of 
precession require uncomfortably large v s , re f, i.e., nowhere close to breakup spin of ~ 1.5 kHz 
for our fiducial mass and radius. (See Cook et al. 1994 for an exhaustive summary of breakup 
frequencies for various equations of state.) Second, if the crust of the neutron star is relaxed 
at its current spin frequency, then P s = l/u S;re f and P p oc v^ef- This re l & tion is shown by 
the solid line in Figure 8. Since we presume that the NS had a relaxed crust in the past, and 
has since then spun down, theoretically feasible precession candidates must lie to the right 
of the solid line. All precession candidates except for the Crab pulsar satisfy this constraint. 
The Crab's claimed precession frequency can also be explained in terms of its crustal rigidity, 
if one assumes that its mass and radius are somewhat different from our fiducial values (e.g., 
M = 1.8M and R = 8.0 km). 
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8. Summary 

A relaxed, self-gravitating object with a rigid crust has a portion of its spin-induced 
bulge that cannot follow changes in the direction of the solid's rotation vector. Such an 
object can precess, at a frequency that is determined by the material properties of the crust 
and the strength of gravity, whose relative strengths are determined by the rigidity parameter 
b. In this paper showed that a NS's rigidity parameter b is a factor ~ 40 times smaller than 
the estimate of Baym & Pines (1971; see Eq. 3). Using this result, we showed that the 
precession frequency v v for PSR 1828-11 is ~ 250(511d/P p ) times faster than expected for a 
NS with a relaxed crust, indicating that its crust is significantly strained. The strain would 
naturally arise from the secular spin down of the star. Large quakes that partially relax the 
accumulated strain energy could also excite precession. Applying our results to PSR 1828-11 
implies the average value of spindown strain in the crust is is a*% e ~ 5 x 10 -5 if the core 
and crust are effectively decoupled, and larger by a factor of ~ 100 if the crust and core 
are dynamically coupled over a spin period. It is not unreasonable to expect the NS crust 
to be able to sustain such modest strain values. Applying our model of a relaxed crust to 
other precession candidates, we found that only SN1987A and the Crab pulsar are consistent 
with the hypothesis of a relaxed crust. The other candidates can be reasonably explained as 
having strained crusts, though unfortunately this explanation has no predictive power; the 
one free parameter (the reference spin) is adjusted to fit the precession timescale. 

We conclude here by mentioning one other application of our result for b. Cutler (2002) 
has shown that for rapidly rotating NS's with relaxed crusts and a strong, interior toroidal 
magnetic field B t , the prolate distortion of the star induced by B t dominates over the oblate- 
ness frozen into the crust for B t > 3.4 x 10 12 G(z/ s /300Hz) 2 . In this case, dissipation tends 
to drive the magnetic symmetry axis orthogonal to the spin direction, and the NS becomes 
a potent gravitational wave emitter. Had one used the Baym-Pines value for b, one would 
have arrived at a substantially larger requirement on the toroidal field (in order for it to 
dominate over crustal rigidity): B t > 1.4 x 10 14 G(z/ s /300Hz) 2 . 

C.C.'s work was supported in part by NASA Grant NAG5-4093, B.L.'s work was sup- 
ported in part by NSF Grant AST 00-98728, and G.U.'s work was supported by a Lee A. 
DuBridge postdoctoral fellowship at Caltech. 
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A. Approximate b from integral of crustal stresses 

Our solution for the NS's rigidity parameter b in § 3-4 involved calculating the exterior 
I = 2 piece of the gravitational potential 5<p for two stars spun up to angular value Q: i) a 
star that is completely fluid, and, ii) a star with relaxed, spherical crust that is otherwise 
identical to the first one. The relative difference in the exterior 5<p is b. Our result (Table 
1) is that b ~ 2 x 1CT 7 , so an accurate calculation of b by the above subtraction method 
requires solving for 5<p to ~ 9 decimal places. Here, as an additional check on our work, we 
estimate b in a way that does not require such high accuracy from our solutions. 

UCB derived an identity relating a NS's mass multipole moment Q 22 to an integral of 
stresses in the crust. The same calculation can be repeated for the Q20 multipole moment. 
Defining Q 2 o by 

Q 20 = J 5p(r, 9, (f))Y 20 (e, 4>)dV (Al) 

we obtain 



Q 20 = -(1-F)- 1 / — [ I ( 4 - ^) *rr + I (6 - C^) *A 



(A2) 



+ 




1 ~. 



1 dU 



dr, 



where g(r) = Gm(r)/r 2 , U(r) = 47rpr 3 /m(r) (note U(r) << 1 in the crust), and the stress 
components of the stress tensor t a b are defined by 



tab = t rr (r)Yi m (r a f b - -e ah ) 

1 

+ t r _i_(r)f ab + t A (r)(A a5 + -e a6 ), 



(A3) 



where the tensors e a b, fab, and A ab were defined in Eq. (59). 

The factor (1 — F)^ 1 in Eq. (A2) requires some explanation. UCB derived Eq. (A2), 
without that factor, within the Cowling approximation, in which the self-gravity of the 
perturbation is neglected. So (1 — F)^ 1 is a correction factor that accounts the effect of the 
perturbation's self-gravity. UCB derived an approximate expression for F (UCB, Eq. 72), 
leading to the estimate that F ~ 0.2 — 0.5, depending on the particular NS model. 

We have calculated the right-hand side of Eq. (A2) for the stresses obtained by solving 
our problem (ii), with the source term (RVt/c) 2 set equal to —1. The result is the residual 
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Q20 for a NS whose crust was relaxed at spin Q = c/R, and which was then spun down 
to v s = 0. We find Q 2 oUiduai = (1 - F) _1 4.24 x 10 38 g-cm 2 for our fiducial 1.4M NS, 
constructed with the AU EOS. 

The rigidity parameter b is the residual Q20 of the non-rotating NS, divided by the Q20 
of the rotating model. The latter is given by 

n 5 ° 2R3 1 / • . % 

V20|rotating — ^ Q Z 5\r=R, (A4J 

where z 5 = 5(p/c 2 . For Q set equal to c/R, we find Q20 1 rotating — 2.8 x 10 45 g-cm 2 for the same 
fiducial model. Thus integrating the stresses in the crust yields b = (1 — F) _1 x 1.5 x 1CT 7 , 
or, for the F in the range 0.2 — 0.5, b ~ 1.9 — 3.0 x 10~ 7 . This is in excellent agreement with 
the value b = 2.47 x 10~ 7 in Table 1. 



B. Analytic solution for b for homogeneous sphere with thin crust 

A classic problem solved by Lord Kelvin was the determination of b = AI^/AIq for the 
case of a constant-p, constant-/^ sphere; the result can be written as b = ^/iV/\E g \, where V 
is the star's volume and E g is its binding energy. (For the derivation, see the classic treatise 
by Love 1944.) In Fig. 9 we plot d E strain /dr for this case. We see that the strain energy is 
concentrated toward the center of the star, which already suggests that this model problem 
will not be a reliable guide for estimating the rigidity of a realistic NS (where all the shear 
stresses are near the surface.) 

Here we solve the corresponding problem for a star with a thin crust. That is, we consider 
a constant-p star that consists of two pieces: a fluid interior (where the shear modulus p is 
zero), plus a thin crust where p is constant. Our case is clearly much closer to a realistic NS, 
where the ratio (crust thickness)/(NS radius) is ~ 1/20. Most of the analysis we require has 
already been carried out by Franco, Link & Epstein (2001), who solved for the displacements 
and strain build-up that occur when such a star spins down. Franco, Link & Epstein (2001) 
work within the Cowling approximation-i.e., they neglect the gravitational perturbation 8(f) 
induced by crustal distortions-so for convenience in this Appendix we neglect Sep as well. As 
explained in Appendix A, the Cowling approximation underestimates b by a factor m 5/3. 

Franco, Link & Epstein (2001) did not restrict themselves to the case of a thin crust, 
but we do so, for convenience. Rather than repeat their derivation, we shall simply quote 
their formulae; e.g., we shall refer to their 39th numbered equation as (F39). Franco et al. 
define R' and R to be the radii at the bottom and top of the crust, respectively. We define 
AR = R — R! and neglect terms that are of quadratic or higher order in AR. We just 
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highlight the basic steps and refer the reader to Franco et al.(2001) for further details. 
We will show that for the thin-crust case, and within the Cowling approximation, 

, 12 fiV c 



(Bl) 



where V c is the volume of the crust and E g is still the binding energy of the entire star. 

When the star's spin-squared changes by 5(£l 2 ), the crust undergoes a displacement 
u(r,9), whose radial piece can be shown to have the form u r (r,9) = f(r)P 2 (9), where P 2 is 
the second Legendre polynomial, and where f(r — R) has the form: 

f(R) = 5 -R 6 -^±(l-b). (B2) 
D v k 

Here v\ = GM/R and M is the star's mass. Eq. (F31) shows that the Eulerian change in 
the (1 = 2 piece of the) gravitational potential exterior to the star is oc f(R), which makes 
it clear that the b on the right-hand side of Eq. (B2) must also equal AI^/AIq. Expanding 
(F23) and (F24) to linear order in AR and subsituting into (F19), we find 

f(R) = l -AR* - ^BR- 2 (B3) 
f'(R) = IaR 2 + ^BR- 3 (B4) 
f'(R - AR) = f'(R) - AR[-2AR - ^BR~ 4 ] (B5) 

where /' = df jdr and A and B are coefficients (to be solved for) appearing in the expansion 
of f(r); see F19. Expanding F38 to linear order in AR and comparing with Eq. (B5), we 
can solve for B in terms of A: 

BR' 3 = —AR 2 + -AR AR . (B6) 
3 5 

Plugging (B6) into Eqs. (B3)-(B4) yields 

f(R) = ^AR 3 -^AR 2 AR (B7) 
f(R) = - l -AR 2 + ^ARAR, (B8) 
and plugging Eqs. (B7)-(B8)into (F33) yields 

f{R) -6 R -^f {1 -n^ ) - (B9) 
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where v\ = /i/p. The second term in parentheses on the righ-hand side of Eq. (B9) is 6, and 
can be re-expressed as 



Note that the coefficient 12/11 in Eq. (BIO) is smaller than the 57/11 in the Baym-Pines 
estimate, Eq. (3), by 209/40 ~ 5.23. So from our uniform, thin-crust model, we see that the 
Baym-Pines expression significantly overestimates the rigidity of a NS. 

If we were to simply replace the 57/11 in the Baym-Pines estimate by 12/11, we would 
still overestimate the rigidity of a realistic NS by a factor ~ 8. To understand this factor, 
it is instructive to to evaluate Q20 (residual for the uniform, thin-crust case using our integral 
expression Eq. (A2), and compare with the realistic case. For the uniform, thin-crust case, 
U — 3, and using Eq. (76) in UCB we find that the correction factor (1 — F)^ 1 is exactly 
5/2. The results in Franco et al.(2001) and a few pages of algebra suffice to show that, to 
leading order in AR, 



for the case of star spun down from Q to zero angular velocity. Also, for a uniform star 
rotating at angular velocity Q, 



Using these results and Eq. (A2), we find that the term oc a rr in Eq. (A2) contributes 
(—^)^fiV c /\E g \ to b, while the term oc a a contributes (|) j^fjV c /\E g \. 

Comparing these intermediate results to those for the realistic NS case (with compress- 
ible matter and a steep density gradient), we see that in the realistic case, there is much 
more cancellation within the integral for Q20. As shown in Fig. 10, in the realistic case both 
o rr and o"a switch signs at different depths inside the crust (so that the contributions to b 
from different layers tend to cancel each other, unlike in the uniform, thin-crust case). Also, 
the a rr and a\ contributions are clearly much closer in magnitude (though still of opposite 
sign) in the realistic case, and so cancel each other much more nearly. 




(BIO) 




(Bll) 
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Fig. 2. — The ratio of the shear modulus to the pressure as a function of density for the NS 
models used in this paper: model AU (solid line), WS (dotted line), nlpoly (dashed line). 
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Fig. 3. — Crustal displacements for a relaxed, nonrotating NS that is spun up to f2 = c/R. 
The different curves represent different equations of state: model AU (solid lines), WS 
(dotted lines), and nlpoly (dashed lines). For v s ,ref = 40 Hz, the values in this plot must 
be multiplied by 6.5 x fCT 5 . Top panel: I = 2 radial displacement z\\ middle panel: / = 2 
transverse displacement z 3 , bottom panel: I = radial displacement y\. 
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Fig. 4. — Radial dependence of the (spindown portion of the) strain tensor for model AU 
(solid line), WS (dotted line), and nlpoly (dashed line), normalized to Q = c/R. For 
v s ,ref = 40 Hz, the values in this plot must be multiplied by 6.5 x 1CT 5 . Top panel: |a rr |, 
second panel: |ay_i_|, third panel: |<ta|, bottom panel: |o"o|. The signs of the components are 
as follows. At the bottom of the crust, for all 3 EOS: a r j_ > 0, o"a < 0, and a > 0. o rr > 
at the bottom for the AU and WS EOS, but a rr < there for the nlpoly model. The sharp 
dips in the figures are zero-crossings. 





0.5 1 1.5 2 

-log 10 (l-r/R) 

Fig. 5. — The distribution of (spindown portion of the) crustal strain a on a slice in a 
meridional plane, containing the rotation axis (pointing upward) for model AU. Results are 
normalized to Q = c/R. Grayscale indicates log 10 er. 
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Fig. 6. — The distribution of (the spindown portion of the) crustal shear stress f on a slice 
in a meridional plane, containing the rotation axis (pointing upward) for model AU. Results 
are normalized to Q = c/R. Grayscale indicates log 10 f, where f is in cgs units (erg cm" 3 ). 
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Fig. 7. — The ratio of energy in crustal stresses to that in centrifugal stresses, 6, scaled to 
the Baym-Pines estimate, plotted as a function of the location of the crust-core boundary 
for equation of state AU (solid line), WS (dotted line), and nlpoly (dashed line). 



-40- 




-3 -2 -1 I 

10 10 10 10 



Spin Period (seconds) 

Fig. 8. — Summary of precession candidates and limits on their initial spins. Squares mark 
spin and precession periods of candidates discussed in § 6, while triangles with arrows indicate 
lower limits on precession periods for B2217+47 and B0959-54. The dashed lines indicate 
minimum precession periods of neutron stars whose crusts last relaxed when they were 
spinning at 1 ms (lower line) and 20 ms (upper line). For a given initial spin period, all 
precessing neutron stars must lie above the corresponding dashed line. The solid line gives 
the precession period for NS's with relaxed crusts. NS's that are slowing down should lie to 
the right of this line. 
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Fig. 9. — The strain energy in a homogenous, elastic sphere, as a function of r. Arbitrary 
units. 
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Fig. 10. — Top panel: pieces of the integrand of Eq. (A2) proportional to t rr (dotted line), t r ± 
(dashed line), and t\ (dot-dashed line), computed for model AU. Bottom panel: cumulative 
integral. 
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Table 1. Summary of NS models, computed b values, and crustal strains 



EOS 


R 


Io 


Ic/Io 


AIn/Io 


b 


b 




^polc 


Cave 




(km) 


Xl0 45 g cm 2 


xlO" 2 


x(y s /kHz) 2 


(actual) 


(B&P) 








AU 


12.13 


1.11 


0.86 


0.259 


2.47xl0~ 7 


9.20x10" 


■6 


17.1 


0.457 


UU 


12.78 


1.23 


1.11 


0.302 


4.20xl0~ 7 


1.25x10" 


5 


18.9 


0.495 


FPS 


11.97 


1.05 


0.78 


0.241 


1.70X10- 7 


8.09x10" 


■6 


16.7 


0.432 


ws 


12.05 


1.10 


0.70 


0.257 


1.47xl0~ 7 


7.62x10" 


■6 


14.9 


0.449 


nlpoly 


12.53 


1.15 


2.86 


0.275 


1.31X10" 5 


2.14x10" 


■4 


103 


0.502 



Note. — Summarizes our results on the crustal rigity and spindown strain for 5 different EOS, and compares 
to the Baym-Pines estimate for b. The average value of the spindown strain <r£„ e and its maximum value o-^ ;e 
(attained on the spin axis, near the North pole) are both normalized to RQ/c = 1. Results here assume the 
density at the bottom is pt, = 1.5 x 10 14 g/cm 3 . 
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Table 2. Precession candidates: observations and inferences 



Name 


Ps 


Ps 








s,ref 


Crustal Strain 




(<0 






(days) 


(Hz) 


(ms) 


a ad , 

pole 


u ave 


B1828-11 


0.405 


6xl0- 14 




511 


40 


26 


2xl0 -3 


5xl0~ 5 


B1642-03 


0.388 


1.78x10" 


■15 


1250 


25 


41 


7xl0~ 4 


2xl0~ 5 


B0833-45 (Vela) 


0.0893 


1.25x10" 


■13 


165 


32 


31 


10~ 3 


3xl0~ 5 


B0531+21 (Crab) 


0.0334 


4.21x10" 


13 


600 


10 


97 


~ 


~ 


B2217+47 


0.5384 


2.77x10" 


■15 


>7300 


<12 


>84 


<2xl0~ 4 


<4xl0~ 6 


B0959-54 


1.437 


5.14x10" 


■14 


>2500 


<33 


>30 


<io- 4 


<3xl0~ 5 


SN1987A 


2.14xl0~ 3 






0.0116 


597 


1.7 


~0 


~0 



Note. — Summarizes our results for 6 precession candidates. JA, ire / is the crust's reference spin value, 
assuming crustal rigidity sets the observed precession period. The strain values are the current spindown strains 
in the NS crust; i.e. the strains induced when a relaxed star spins down from the reference spin u a ^ re f to the 
current spin v a . For B2217+47 and B0959-54, only lower limits on the precession period are claimed in the 
literature, yielding lower limits on v s ^ re f. For the Crab and SN1987A, u s ^ re f is less than, or comparable to, v s 
for our fiducial NS model, so there is no evidence of nonzero crustal spindown strain. 



